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ABSTRACT 



CoRoT-3b is a 22 Jupiter-mass massive-planet/brown-dwarf object, orbiting an F3-star with a period of 4.3 days. We analyzed the 
out-of-transit CoRoT-3 red-channel lightcurve obtained by the CoRoT mission and detected the ellipsoidal modulation, with half the 
orbital period and an amplitude of 59 ± 9 ppm (parts per million), and the relativistic beaming effect, with the orbital period and an 
amplitude of 27 ± 9 ppm. Phases and amplitudes of both modulations are consistent with our theoretical approximation. 
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1. Introduction 

Close binary stellar systems display two well-known periodic 
photometric modulations — the ellipsoidal effect, due to the dis- 
tortion of e ach compone nt by the gravity of its companion (see 
a review by|Mazeh 2008), and the reflection/heating effect (re- 
ferred to here as the reflection effect), induced by the luminosity 
of each s tar, which affec ts only the close side of its compan- 
ion (e.g., iFor et alj|2010b . These two effects can be observed 
even for non-eclipsing binaries, but are much easier to study 
in eclipsing binaries, where the binarity of the system and the 
phases of the orbital motion are well known from the observa- 
tions of the eclipses. Most algorith ms that analyze lightcurves of 
eclips ing binaries, such as EBOP (Etzel 1980; Poppe rand Etzell 
119811) and its derivative EBAS (Tamuz, M azeh & North 2008), 
WD (IWilson & Devinnevlll97lh . and ELC (Orosz & Hauschildt 
2000), include by default these two effects in their model of the 
out-of-eclipse lightcurve. 

A much smaller and less studied photometric modulation is 
the relativistic beaming effect, sometimes also called Doppler 
boosting, induced by the stellar motion relative to the observer 
— Vreh whose amplitude is on the order of V re i/c, where c is the 
velocity of light. Before the era of space photometry this effect 
has been noticed only once, by Maxted, March & North (2000), 
who observed KPD 1930+2752, a binary with a very short pe- 
riod, of little longer than 2 hours, and a radial-velocity amplitude 
of 350 km/s. The beaming effect of that system, which should be 
on the order of 10~ 3 , was hardly seen in the photometric data. 

Space photometry, which was developed to detect the minute 
transits of exoplanets, has substantially improve d the preci- 
sion o f the produced lightcurves. The CoRoT dRouan et al.l 
Il998t iBadin et alJ l2006t lAuvergne et all 120091) and Kepler 
dBorucki et al.ll2.010l iKoch et al.ll2010h missions are producing 
hundreds of thousands of continuous photometric lightcurves 
with timespan of tens and hundreds of days, with precision 
that can reach as high as 10~ 3 -10~ 4 per measurement. It was 
therefore anticipate d that CoRoT and Kepler should detect all 
three effects (e.g., [Drake 2003), in particular the beaming ef- 



fect for both planets (Loeb & Gaudi 2003) and eclipsing binaries 
(Zucker, Mazeh & Alexander 2007). 

As predicted, Ivan Kerkwiik et all d2010h detected in the 
Kepler lightcurve the ellipsoidal and the beaming effect of two 
eclipsing binaries, KOI 74 and KOI 81 (Rowe 2010). They used 
the radial-velocity photometric beaming effect to derive the mass 
of the secondary in the two system s and showed that in both 
cases it was probably a white dwarf. IWelsh et"aD (l2010h identi- 
fied the ellipsoidal effect in the Kepler data of HAT-P-7, a system 
with a kno wn planet of 1 .8 J upiter masses (=M./ t lp ) and a period 
of 2.2 days dPal et al.ll2008l) . ISnellen et al. (2009) detected in the 
CoRoT data the reflection effect of CoRoT- 1 . 

In this paper, we report the detection of the ellipsoidal 
and the beaming effects of CoRoT-3, indu ced by its massive- 
planet/brown-dwarf companion. CoRoT-3b dDeleuil et al.| [2008) 
is a 22 Jupiter-mass object, orbiting an F3-star with a period 
of 4.26 days. The stellar rotation is probably synchronized with 
the orbital period. We analyzed the CoRoT-3 out-of-transit red- 
channel lightcurve and detected two modulations, one with the 
orbital period and the other with its first harmonic. We attributed 
the two modulations to the beaming and the ellipsoidal effects, 
respectively, as their phases and amplitudes were consistent with 
our order-of-magnitude approximation. Sect. 2 presents our data 
analysis, Sect. 3 compares our findings with theoretical approx- 
imations, and Sect. 4 summarizes our results. 



2. Data analysis 



CoRoT-3 was discovered (Deleu iTet alJ l2008h in the data ob- 
tained during the first long run of the CoRoT mission — LRcOl, 
which lasted for 152.012 d, fro m May 26 until Octo ber 25, 2007 
(for details about this run see Cabrer aet al.l 120091) . The optics 
of the mission include a bi-prism that disperses the stellar light 
into three channels, red, green, and blue, the sum of which is 
called the white channel. For bright stars, including CoRoT-3, 
the light intensity coming through each of the t hree channels 
is ava ilable. We used the so-called N2 data level (Baudi n et al.l 
2006) of CoRoT-3, which is now public. 
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As the modulations we searched for were quite small, we 
had to prepare and clean the data before searching for any peri- 
odic effect. We decided to concentrate on the red-channel data, 
because this channel included most of the stellar light detected 
by CoRoT, and th e other two c h annel s only added noise to the 
data (see below). ISnellen et al.l (j2009) adopted a similar strat- 
egy when analyzing the lightcurve of CoRoT-1. This section de- 
scribes how we 'cleaned' the data, removed the long-term vari- 
ation, and searched for the periodic modulation with the orbital 
period and its harmonics. 

2.1. Cleaning the lightcurve 

The 'cleaning' of the CoRoT-3 lightcurve had the following 

stages: 

- Rebinning: Corot-3b was detected before the CoRoT run was 
completed and therefore the cadence of the observations was 
changed during the run — the first part of the lightcurve is 
composed of 5 12 s exposures, while the later part contains 32 
s exposures. Since we were interested in modulations with 
periods equal to or longer than half the orbital period, at 
about 2.1 d, the entire light curve was re-binned into 512 
s bins. Altogether, we derived 22,072 valid measurements. 

- Removing transits: 1165 measurements taken during the 
transits of CoRoT-3 were removed from the analysis. 

- Jump removal: One 'jump' , at CoRoT HJD of 2746.99, prob- 
ably caused by a 'hot pixel' event, was identified and cor- 
rected. The counts after the jump were adopted to the stellar 
flux before the jump, while 24 measurements following the 
jump were removed. 

- Outlier removal: We identified 106 outliers by calculating the 
running median and RMS around each point, and rejecting 
measurements that differed by 4<x or more from their corre- 
sponding median. We were left with 20,801 data points. 

2.2. Long-term detrending with a cosine filter 

The CoRoT-3 lightcurve clearly contained a long-term variation, 
as can be seen in Fig. 1, where we plot the relative red-channel 
flux, after subtracting and dividing the original flux by its me- 
dian. 

To remove this trend we used a discrete cosine transform 
dAhmed et al.l 11974b . adopted to the unevenly spaced data we 
had in hand. We fitted the data with a linear combination of the 
first N low-frequency cosine functions 



f 2n ) I 2T 
fj(tj) = cos( — i x tj); i = 0, N > , where N = Round 

\ 2T J \APorb 



T = 152.012 d is the timespan of the observations, P or b = 
4.2568 d is the orbital period of CoRoT-3, and tj is the timing 
of the j-th measurement. The fitting finds the linear coefficient 
di for each of the cosine functions, so that the fitted model is 



M(tj) = Yj atfiitj) 



(2) 



We then subtracted the model MUfj) from the lightcurve. 

The general idea was to perform a high-pass filter, so 
we removed all the low-frequency cosine components of the 
lightcurve without altering the periodic modulation of the or- 
bital period. In a similar manner we also removed the satel- 
lite and earth freq uencies, which appeared in the N2 data (e.g., 
iMazehet al.lf2009h . 



Fig. 1 shows the red-channel lightcurve before and after the 
removal of the long-term trend and the satellite and earth modu- 
lation. The RMS of the cleaned lightcurve is 904 ppm (parts per 
million). A similar analysis of the blue- and green-channel data 
yielded lightcurves with highly correlated noise, of an RMS of 
1600 and 2000 ppm, respectively. The white-channel lightcurve, 
which includes the blue and the green data, is affected by similar 
problems. These results supported our decision to consider the 
red-channel data only. 

2.3. Fitting the amplitudes of the ellipsoidal, beaming, and 
reflection effects 

We finally proceeded to fit a model that includes the ellipsoidal, 
beaming and reflection effects. We approximated each of the 
three effects using pure sine/cosine functions, relative to the mid- 
dle of the transit, t tmn , denoted as phase zero. The reflection and 
the beaming effects were approximated by sine and cosine func- 
tions, respectively, with the orbital period, and the ellipsoidal ef- 
fect by a cosine function with half the orbital period (see next 
section). In this approximation, we expressed the stellar flux 
modulation AF as a fraction of the averaged flux F, and a func- 
tion of t = t - tfron '■ 



AFemp(f) 



= -Aellip COS 



In 



Porb/2 



AF be!im (t) . I 2n 

7; — n-beam Sin — t 

F \P rb 



^ = - Are/ ' C0S fc f 



(3) 



(4) 



(5) 



where the coefficients, A e m p , A^aming, and A re fi are all positive. 

We therefore fitted the cleaned, detrended lightcurve of 
CoRoT-3 with a 5-parameter model, M e b r , consisting of two fre- 
quencies 



In 



2n 



Mebritj) = flo+flic cos \j r ~ tjj+du sin tjj+a 2c cos y p J 



In 



tj \+a 2s 



as performed by ISirko & Paczvnskil ((2003). The fitting process 
could find any value, positive or negative, for the five parame- 
ters. However, we did expect a\ s to represent the beaming effect 
and therefore be positive, a\ c to represent the reflection effect 
and therefore be negative, a%c to represent the ellipsoidal effect 
j arid Aerefore be negative, and a 2s to be close to zero. The last 
expectation is important because it ensures that the ellipsoidal 
modulation is detected with the correct phase. The ao parame- 
ter was induced to remove any DC component left in the data. 
Obviously, the expected absolute values of the three coefficients 
depended on the parameters of the CoRoT-3 system. Thus, after 
the analysis was performed we were able to verify our results by 
comparing them with the expected values. 

The results of the fitting are given in Table Q] (with some 
order-of-magnitude theoretical expectations; see next section) 
and plotted in Fig. [2] In the figure one can easily discern the el- 
lipsoidal modulation, with half the orbital period, and the beam- 
ing effect, which causes the difference between the two peaks. 
One indication of the consistency of our results with the ex- 
pected modulations is the correct sign of the three first coef- 
ficients and that the fourth coefficient is smaller than the third 
one by at least one order of magnitude. This suggests a close 
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CoRoT-3 red-channel lightcurve 




Time (days) 

Fig. 1. The lightcurve of CoRoT-3, before (blue) and after (red) detrending. The long-term model is presented by the wide line. 
Time is the CoRoT time = BJD - 2451545.0. Flux is relative to the mean red-channel flux. The detrended lightcurve is shifted by 
0.02. 



Table 1. The fitted coefficients and the theoretical expected am- 3.1. Ellipsoidal effect 



plitudes of the three effects of CoRoT-3 
Coefficient Derived value Expected amplitude Effect 



(ppm) 



(ppm) 



a 2c 

rtlj 



-14 ±9 
27 ±9 

-59 + 9 
0.1 ±9 



-a nfl g.2 ± 0.3) 
«w(29 ± 0.5) 
-a e m p (32 ± 5) 



Reflection 

Beaming 

Ellipsoidal 



agreement between the orbital phase of the ellipsoidal modula- 
tion and that of the transit. We note that the errors in the four 
coefficients are all 9 ppm, or somewhat smaller than 10~ 5 in rel- 
ative flux. This precision does not allow a significant detection 
of the reflection modulation. However, the detection of both the 
ellipsoidal and the beaming effect is highly significant. A boot- 
strap test indicates that the probability of detecting the beaming 
modulation by chance is 2 X 10~ 4 . 



To derive the order-of-magnitude o f the ellipsoidal effect w e 
used the analy tical appro ximation of Morris & Naftilan| (1 19931) . 
who used the iKonall (119591) expansion of the periodic variation 
into discrete Fourier series with terms that depend on the ratio 
R*/a, where R, is th e stellar radius an d a is the semi-major axis 
of the orbit (see also Pfahl et al. 2008). Assuming R*/ais small, 
the leading term of the stellar variation has a semi-amplitude of 



where 



= 0.15 



(15 + «)(!+*) 



(7) 



(8) 



is of order unity. In the above expression, m p is the planetary 
mass, M» is the stellar mass, g is the stellar gravity da rkening 
coefficient, and u is its limb-darkening coefficient (e.g., Mazeh 
2008). 



3. Theoretical approximation 

This section presents our theoretical approximations of the el- 
lipsoidal, beaming, and the reflection effects. We are not inter- 
ested in detailed calculations, which depend on specific models. 
For example, to calculate the ellipsoidal modulation one could 
integrate the light originating from the individual surface ele- 
ments of the rotating star, which is slightly d eformed by the tidal 
force induced by its small companion (e.g. JOrosz & Haus childt 
2000). Such an analysis was carefully perfo rmed very recentl y 
to model the Kepler lightcurve of HAT-P-7 dWelsh et al .1120101) . 
Instead, we are interested here only in order-of-magnitude ap- 
proximation, so we can check the consistency of the detected 
amplitudes with the theory. All our approximations are evaluated 
for an inclination angle i close to 90 deg, so we could ignore the 
sin i dependence of the three effects. 



3.2. Beaming effect 

For a circular orbit, the amplitude of the beaming effect can be 
written as 

K r 

Abeam — ^beam 4 , (9) 
C 

where K r is th e stellar radial-veloci t y amplitude and c is the 
speed of light dLoeb & Gaudil 120031: IZucker et al] 120071) , and 
<Xbeam is of order unity. The factor 4-K, /c represents the beaming 
effect for bolometric photometric observations, but ignores the 
Doppler shift photometric effect, which appears when the pho- 
tometric observations are made in a specific bandpass, so that 
some of the stellar light is shifted out of or into the observed 
bandpass. The latter is accounted for by the abeam factor, and we 
assume that for the CoRoT red bandpass it is of order unity. 
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Fig. 2. The folded cleaned lightcurve of CoRoT-3, binned into 100 bins, with the fitted model. The residuals are plotted at the bottom 
of the figure 



3.3. Reflection effect 

In our simplistic approximation we include in the reflection 
modulation the thermal emission from the dayside of CoRoT- 
3b, assuming both are modulated with the same phase (e.g., 
ISnellen et al.l 2009). The amplitude of the modulation of the re- 
flected light alone is 



A re fl — Pgeo ' 



(10) 



where r p is the planetary radius a nd p Qeo is the geometri- 
cal albedo (e.g., iRowe et alj |2008|) . iRoweet alj d2008l) found 
quite a small albedo, o f 0.03, for HD 209458, but recent study 
dCowan & Agoll l2009f) suggested that exoplanets may have a 
much larger albedo, of up to 0.5. We therefore write the am- 
plitude of the reflection effect, including the thermal emission, 
as 



i-refl 



trefl 



0.1 



(ID 



where a re fi is of order unity. 
3.4. Application to CoRoT-3 

Table |2] presents the relevant par ameters of CoRoT- 3. The first 
four parameters were derived bv lDeleuil et alj (12008), while the 
last parameter, the radial-velocity amplitude, was deduced by 
iTriaud et al.l d2009l) . From these parameters we derived the ex- 
pected values of the amplitudes of the ellipsoidal, beaming, and 
reflection effects, which are given in Table Q] 

We emphasize that in each of the three theoretical amplitudes 
given in Table [T| the main source of uncertainty is hidden in 
the a factor, which we did not calculate. The numerical values, 
with their relatively small errors, are only order-of-magnitude 
approximations for inclination angles close to 90 deg. 

The amplitudes derived from the cleaned lightcurve of 
CoRoT-3, as shown in Table Q] are of the same order of mag- 
nitude as the expected values, based on our simplistic approxi- 
mation. This is true in particular for the beaming effect, where 



Ta ble 2. CoRoT-3 param eters, as derived bv lDeleuil etall (120081) 
and lTriaud etall (120091) 



Parameter Derived value 



Unit 



a/R, 


7.8 ±0.4 




r p /R, 


(663 ± 9) x Kr 4 




m p 


21.7 ± 1.0 


M Jup 


M„ 


1.37+0.09 


M Q 


T, 


6740 ± 140 K 


deg 


K r 


2170 ±30 


m s _1 



the theoretical approximation was found to be quite accurate. 
We therefore propose that we have detected the ellipsoidal and 
beaming effects of CoRoT-3. Our results suggest that the a e m p 
factor in CoRoT-3 is on the order of 2. The reflection effect was 
too small to ensure a significant detection, given the SNR of the 
lightcurve. 

4. Discussion 

Our analysis has demonstrated that the red-channel lightcurve of 
CoRoT-3 includes the ellipsoidal and beaming effects. This is the 
first time that the beaming effect has been detected for substellar 
companion. We have been able to detect the two effects, with 59 
and 27 ppm amplitudes, respectively, because of a combination 
of three features: 

- CoRoT-3 brightness: With r'-mag of 13.1, the star is among 
the brightest Co RoT targets, which are typically in the range 
of 1 1 to 16 in t> (iDeleuil et al.ll2008l) . 

- The long observational run: The LRcOl lasted for 152 days, 
and the CCDs did not show yet any aging signals. 

- The mass of CoRoT-3b: This massive-planet/brown-dwarf 
companion has one of the largest masses, 22 Mj up , discov- 
ered by CoRoT for substellar objects. 

The last feature suggests that the stellar rotation has achieved 
synchronization with the orbital period of 4.3 days, without 
which the analysis of the ellipsoidal effect could have been more 
complicated. 
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Had this analysis been performed immediately after the dis- 
covery of the transits of CoRoT-3, and in particular before the 
radial-velocity confirmation of the planetary nature of the tran- 
siting object, we could have estimated the mass of the unseen 
object from the observed amplitudes of the ellipsoidal and beam- 
ing effects, provided we had been able to accurately derive their 
expected values. Such analysis could, in principle, save costly 
radial-velocity observations, or at least reduce to a minimum the 
number of observed velocities needed to confirm the substellar 
mass of the transiting object. 

Obviously, the analysis of lightcurves obtained by space mis- 
sions is dramatically different from those obtained by ground- 
based photometry. For the latter, the appearance of the ellip- 
soidal modulation in the data of transit candidates was consid- 
ered a sig n that the transiting o bject was of stellar nature, as sug- 
ge sted by [ Sirko & Paczvhski (20 03b, and applied, f or example, 
bv lKane et alj d2008l) and lPietrukowicz et alj (120101) . This is so 
because of the relatively high threshold of detection of the ellip- 
soidal modulation in the ground-based photometry. On the other 
hand, the present work, and the study of Wel sh et al.l 12010). sug- 
gest that the detection of the ellipsoidal modulation with a small 
amplitude in the CoRoT and Kepler data may indicate that the 
transiting object is a massive-planet/brown-dwarf object. 

The present analysis suggests that, in principle, the three ef- 
fects, or at least two of them, can be detected in the CoRoT 
lightcurves for some massive-planet/b rown-dwarf objects, even 
without any t ransits , as suggested by lLoeb & Gaudil d2003l) and 
IZucker et alJ(l2007h . The effects can be stronger for systems with 
shorter orbital periods, and therefore can be detected in stars 
fainter than CoRoT-3 in the CoRoT fields. Many objects simi- 
lar to CoRoT-3 should also be detected by Kepler, because of 
both higher SNR and longer timespan of its lightcurves. 
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